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ABSTRACT 



By observing radiation- affected gas in the Cepheus B molecular cloud we probe whether the sequential star formation in this source 
is triggered by the radiation from newly formed stars. We used the dual band receiver GREAT onboard SOFIA to map [C ii] and CO 
13-12 and 11-10 in Cep B and compared the spatial distribution and the spectral profiles with complementary ground-based data of 
low-/ transitions of CO isotopes, atomic carbon, and the radio continuum. The interaction of the radiation from the neighboring OB 
association creates a large photon-dominated region (PDR) at the surface of the molecular cloud traced through the photoevaporation 
of C^. Bright internal PDRs of hot gas are created around the embedded young stars, where we detect evidence of the compression 
of material and local velocity changes; however, on the global scale we find no indications that the dense molecular material is 
dynamically affected. 
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1. Introduction 

The molecular cloud Cepheus B, located at a distance of about 
730 pc can be regarded as a prototype of a region with sequen- 
tial star formation. It has a wedge-like shape, pointing NW to 
the young Cepheus OB 3 association, including the luminous 
stars HD 217086 (07) and HD 217061 (Bl). The OB associ- 
ation is composed of two subgroups of diff'erent age s, where 
the youngest group lies closer to the molecular cloud ^Sargent' 
|1977| ). The radiation from the association creates a bright H ii re- 
gion, S 155, at the cloud s urface. Deeper in t he cloud, star forma- 
tion contin ues (Moreno-Corral et al.||1993| ), creating a hot core 
( |Felli et aL] 1978 ) that was resolved int o two ultracompac t Hii 
regions and one extended Hii region by Test! et al. ( 1995| ). The 
main energy source in the hot core is a B 1 star in source A of 



[Test! et aL] ( |1995 ), but m ore pre-main- sequence stars are con- 
tained in the hot core (Moreno-Corral et al. 1993| ). The density of 
molecular material increases further towards the SE, where two 
neighboring column density peaks have bee n identified through 
observations of low-/ ^^CO and C^^O (Ung erechts et al.||2Q00 



Mookerjea et al.||2QQ6| ). The molecular cloud observations by 
Beuther et al.H2000| ) showed a clumpy structure that allows a 
deep penetration of UV radiation from the OB association and a 
large-scale velocity gradient with blueshifted velocities towards 
the OB association. At a resolution of 43'', three [Cii] s pectra at 
158 jim we re observed in Cepheus B using the KAO ( [Boreiko 



et al. II 1990 1. The [Cii] spectra was found to have two velocity 



components, one associated with the molecular material and one 



with S I 55. The CO 4-3 and [C i] observations of Moo kerjea etal 
( 2QQ6| ) showed that the main heating of the molecular material 



occurs in PDRs (photon-dominated regions) created below S 155 
and around the hot-spot H ii regions. The FUV field is estimated 



to be between 1500 to 2 5;^o, where ;^o is the average interstellar 
radiation field ( |Draine|1 978). 

To look for indications of sequential star formation due to ra- 
diative triggering as proposed by |Getman et al.|p009 l, the den- 
sity, temperature structure, and kinematics of the FUV irradiated 
PDRs at the cloud surfaces need to be probed. Signs of trigger- 
ing processes would be the formation of density enhancements 
correlated with the radiative pressure increase at the cloud sur- 
face (radiative shells, see e.g. Dale et al.|2009| ) or indications of 
velocity gradients in the molecular material created below the 
radiative compression zones. 

To address these questions, we performed spectral mapping 
observations of Cep B in [Cii] and CO 13-12 and 11-10 using 
the Germ an REceiver for Astro nomy at Terahertz Frequencies 
(GREAT i Heyminck et al.| [2012f] onboard the Stratospheric 
Observatory For Infrared Astronomy (SOFIA). The new data 
were combined with existing observations of low-/ CO iso- 
topologs and [C i] to probe the structure of the PDR and the 
velocity distribution. 



2. Observations 

We observed the hot core region of Cepheus B during the nights 
of April 7 and July 19, 2011 with GREAT on SOFIA at an alti- 
tude of ~ 42000 feet. During the two nights were observed ar- 
eas of the same size (2'4xlf6), but off'set relative to each other. 
GREAT is a modular heterodyne instrument, with two chan- 
nels used simultaneously. In April, channel Ll#a was tuned to 



^ GREAT is a development by the MPI fiir Radioastronomie 
and KOSMA/Universitat zu Koln, in cooperation with the MPI fiir 
Sonnensystemforschung and the DLR Institut fiir Planetenforschung. 
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CO 11-10 (1267.015 GHz) and channel L2 was tuned to [Cii] 
^P3/2 ^^Pi/2 (1900.537 GHz) for four coverages of the region. 
In July, Ll#b was tuned to CO 13-12 (1496.923 GHz), while L2 
was tuned to [Cii] for three coverages. The HPBWs are 2 IV 3 at 
1.3 THz, 19^.^6 at 1.5 THz, and 15'' at 1.9 THz. 

The observations were performed in position- switched on- 
the-fly (OTF) mode with an OFF position at ^a = -600''; AS = 
-h600''. As backend we used the acousto optic spectrometer with 
a bandwidth of 1 GHz and a resolution of 1.6 MHz. The integra- 
tion time was 0.6 minutes/position. The measured system tem- 
peratures were 3100K for CO 11-10, 2500K for CO 13-12, and 
3200-3500 K for [Cii]. The forward efficiency is T]f = 0.95 and 
the main beam efficiencies (^mb) are 0.51 for [Cii] and 0.53 for 
CO 11-10 lines ( |Guan et al. |2012|). Details of calibration of the 
data are presented by Guan et al. p012 ). We fitted and sub- 
tracted baselines of the fourth or sixth order to the [C ii] and CO 
data. 

Complementary observations of [Ci] (^Pi ^^Po at 
492.161GHz) were performed in December 2002 using the 
CHAMP 2x8 pixel array receiver ( [Guesten et al.||1998| ) at the 
CSO. The fully sampled [Ci] map of 3'3x3f7 was observed 
in position-switched mode (reference position -300'', 0). The 
HPBW for [C i] was 14'.'5 at 492 GHz and tjmb was 0.51 ( |Philipp 
et al.||2006 ). The [Ci] data were smoothed to 20" for compari- 
son. 



3. Morphology 

Figure [T] shows the integrated intensity maps of [Cii] and 
CO 11-10 and CO 13-12. The (0,0) position of the maps is at 
a = 22^57'^8'.7 and S = 62°34'23" (J2000). The [Cii] emission 
shows a single peak (-34'.' 7, -22'.' 9) located to the W of an em- 
bedded Bl star, and the emission drops towards the E. The CO 
11-10 exhibits three peaks referred to as #1, #2 and #3 in the rest 
of the paper, forming a shell around the [C ii] peak. The emission 
is extended further in the NW-SE direction. 

Figure [2^ compares the [Cii] and CO 11-10 observations 
with radio continuum observations by| Testi et al .|(|1995|). Testi 
et al. show that the 8.4 GHz radio continuum emission compo- 
nent A, a blister H ii region, peaks slightly to the SW of the posi- 
tion of the B 1 star. We find a sequence of peaks in that direction 
starting with the radio continuum close to the heating star, the 
[Cii] peak, and the CO 11-10 peak #2, consistent with the ex- 
pected layered structure of a PDR. The CO 11-10 peaks #2 and 
#3, both S of the source A trace dense and w arm molecular gas . 
Peak #1 coincides with the radio source B ( [Testi et ar]|1995| ). 
The CO 13-12 emission appears to be centered close to the Bl 
star, but the detailed distribution around it could not be observed 
due to time constraints. 

Figure [2p compares our data w ith IRAM 30m ^ ^CO 2 -1 ob- 
servations at 10" resolution by Un gerechts et al.| ( |2000| ), trac- 
ing molecular gas with high column densities, which has a ridge 
like structure extended towards the main cloud E of the map. 
A similar intensity distribution is also seen in the ^^CO 1-0 map 
( Ungerechts et al.|2000|), the [C i] 1-0 map (see Fig. [3]), and in the 
CO 4-3 map (at 1' [ Mookerjea et al.|2006 l). The CO 11-10 peak 
#3 coincides with a column density peak in these maps, while 
peaks #1 and #2 indicate an increased CO excitation, e.g. due to 
a higher temperature of the molecular material. More PDR lay- 
ers are formed at the surface of the molecular material in the SE. 
They are also traced by the [C ii] and the high-/ CO lines but at 
a low level of emission. We detect almost no CO emission to the 
N and NE of the embedded B 1 star. 
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Fig. 1. (a) Velocity-integrated intensity map of [Cii] in color 
and contors. The contor levels are at 79 to 339.4 Kkms"^ 
(peak) in steps of 29 Kkms"^ (b) Integrated intensity map of 
CO 11-10 (color) overlay ed with contors of CO 13-12 (white) 
and [Cii] (black) intensities. Contour levels for CO 13-12 are: 
2.1 to 11.7 Kkms"^ (peak). VLSR-range of integration is [Cii]: 
[-20; -5], CO: [-16; -11] kms"^ The white star denotes the 
Bl star at (a;6) = (22^57'^6.2^62°37'55.4"). Off'sets of the 
marked CO peaks #1, #2, and #3 are (-45",-h64"), (-41",-h20"), 
and (-9",-f3"), respectively. 




Fig. 2. (a) Optical image (grayscale) overlayed with contors of 
8.4 GHz continuum emission (black) ( [Testi etal.|1995[ ) and con- 
tors of [C ii] (green dashed) and CO 11-10 (red), (b) Intensity of 
^^CO 2-1 emission integrated between vlsr = -17 to -10 km s"^ 
( Ungerechts et al. ''2000) (grayscale) overlayed with contors of 
CO 11-10 (red) and [C ii]_( green dotted). Contour levels of [C ii] 
are the same as in Fig.[l] for CO 11-10: 6.8 to 32.1 Kkms"^ 
(peak) in steps of 3.6 Kkms"^ The star and the triangle denote 
the positions of Testi's NIR sources A and B, respectively. 
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Since [C ii] can also originate in ionized gas, it is reasonable 
to expect enhanced [C ii] emission at the location of the radio 
continuum peaks detected by Testi et al.| ([T995). However, the 
[C ii] emission shows no enhancement at the location of the CO 
peak #1, which coincides with the B -component. This indicates 
that the source creating the B -component of the radio contin- 
uum emission ionizes at most a small fraction of the molecular 
material, consistent with Testi 's suggestion that source B is a 
protostar, not yet ionizing, but heating up its surroundings. 

Our observations are consistent with Testi 's suggestion of a 
blister-type H ii-region for component A, created by the B 1 star, 
and a density gradient of the surrounding material leading to an 
"open cone" configuration. The density decreases to the N, so the 
cone opens in this direction. The ionized material may escape 
to the N, but cannot expand much to the S because of denser 
molecular material there. This results in increased and sharply 
bounded radio continuum emission to the S and in extended, but 
lower intensity emission to the N of the position of the embedded 
star. 



4. PDR layering 

As discussed above we find a layering sequence of [C ii], high-/ 
CO and low-/ CO peaks around the H ii region associated with 
the Bl star. Figure [3] shows the eff'ect of this layering in terms of 
the mutual shift of the peaks of the [Cii], CO 11-10, and ^^CO 
2-1 emission. We find a separation of about 0.02 pc between 
[C II] and CO 1 1-10 and of 0.08 pc between [C ii] and ^^CO 2-1. 
For [C i] the peak emission is seen at the edge of the mapped 
region. The [Ci] distribution is similar to that of ^^CO 2-1, but 
the emission is weak in the considered cut. The shift in the peaks 
is consistent with the predictions of a uniform PDR model, al- 
though in reality, it will also be aff'ected by small-scale density 
variations. 

For a qualitative com parison we used the KOSMA-r PDR 
model (Rolli g et al.|2006| ) to simulate gas temperature and chem- 
ical structure inside the PDR. KOSMA-r considers a spherical 
clump illuminated by an isotropic FUV radiation field and cos- 
mic rays. Input parameters are the radiation field, the clump mass 
and the density at the clump surface. It assumes that the clump 
has a power-law density profile of ^(r) oc r~^-^ for 0.2 < r/rd < 1 
and n(r) = const, for r/rd ^ 0-2. To simulate a plane-parallel 
PDR through a spherical model, we chose a clump of 10^ M©, 
i.e. with a mass much larger than the available molecular mass 
(see [Mookerjea et al.|[2006| ). The outer layers of this massive 
clump then form an almost plane-parallel configuration suitable 
to studying the observed stratification. 

Figure [4] shows the local optically thin emissivity, i.e. the 
intensities per cloud depth without optical depth correction, for 



two models with a density n(rci = 10^ cm ^ and n(rci = lO^cm 



illuminated by a radiation field of 1000;^o (Mook erjea et al. 
|200 6). The model predicts a stratification with a 0.001. . . 0.01 pc 
separation of the CO 11-10 peak relative to the [Cii] peak and 
a further separation of 0.12 . . .0.3 pc to the ^^CO 2-1 pea]^ 
indicating that the parameters used are reasonable, but do not 
provide a full match to the observations. The observed intensity 
ratios of the three lines suggest that the density in the cloud is 
closer to 10^ cm"^ than to 10^ cm"^. Since the CO 11-10 line 
is very sensitive to the gas density, the similar peak intensities 



^ ^^CO 2-1 always shows a broad peak filling a large fraction of the 
clump core, independent of the clump mass, indicating that the displace- 
ment of the ^^CO 2-1 peak relative to the PDR surface is fairly well 
determined by the density structure than by illumination. 



at #2 and #3 indicate approximately equal gas densities, while 
the higher ^^CO and [Ci] intensities at peak #3 indicate higher 
column densities of molecular gas there. 

The model, however, predicts a narrow, but weak [C i] peak 
in front of the CO 11-10 layer, while the [Ci] emission is very 
extended on the map [^ In contradiction to the model, the [C i] 
emission seems to measure the molecular column density, so not 
following the stratification of the other tracers. 
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Fig. 3. Left: Map of integrated intensity distributions of ^^CO 2- 
1 (greyscale, red contor), [Cii] (black), CO 11-10 (blue), and 
[Ci] (green). Right: Cut through the [Cii] peak and the CO 11- 
10 peak #2. The black dotted vertical line denotes the position 
closest to the Bl star along the cut. Integration ranges: [C ii] and 
[Ci] [-20; -5], CO 11-10, and ^^CO [-17; -10] (kms"^). 



CI] 1-0 (xioy ^ 




(a) ; 



1.5 10 

2 10* r 

..5 10* - 

10*^ 

5000 - 

- 

10 



[Cll 1-0 (x50j 
to 11-10 (x500) 
"cO 2-1 



.(b) 




Fig. 4. Optically thin line intensities per cloud depth for the two 
model densities: (a) ^surf = 10^ cm"^, (b) ^surf = 10^ cm"^. 
Radiation field and clump mass are 1000;^o and 1000 M©, re- 
spectively. 



5. Velocity structure 

To study the main characteristics of the [C ii] velocity distribu- 
tion we show in Fig. [5] the map of the centroid velocities and 
the corresponding position-velocity diagram for a cut through 
the map. The map of the centroid velocities was generated by 
fitting a single-component Gaussian profile to each spectrum. 
The velocity distribution is dominated by a large-scale gradient 
of about 1 kms~^ over 200'' from NW to SE, visible as a color 



^ Spot checks with the PDR models by M. Kaufman and the Meudon 
group show the same inversion in the [C i] layering that we find with the 
KOSMA-r model. 
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gradient in Fig.|5^ and as a global slope in Fig. [5}). This is con- 
sistent with the velocity gradient of 4.2 x 10"^ kms"Varcsec in 
the molecula r cloud derived from the channel maps of CO 3-2 
observations ( Beuther et al. 2000 ). Moreover, one can see a good 
correlation between the centroid velocity map and the ^^CO map 
in Fig.[2]in the sense that pixels associated with strong emission 
in low- and mid-/ CO have systematically higher velocities than 
nearby pixels measuring lower density material. The velocity of 
the low-density material is closer to that of S155; i.e., the ion- 
ized gas is always somewhat blueshifted relative to the denser 
molecular cloud. 

The position- velocity diagram Fig. [5})) shows a broadening 
of the [C ii] line at the position of the [C ii] peak and there ap- 
pears to be an additional velocity component visible as a blue 
shoulder. An interesting feature is, however, the change in the 
velocity gradient around the [C ii] peak. The contors in the PV 
diagram form a parallelogram NW of the peak where the large- 
scale velocity gradient is inverted. The inversion is not very pro- 
nounced when considering the Gauss-fit velocities along the cut, 
but is apparent as higher velocity in the map (Fig. [5^) around CO 
peak #2. 

To look at the details of the spectral profiles, we averaged the 
spectra around the peak positions #1, #2, and #3 over an area of 
33'' diameter (Fig.|6]). We find emission from diff'erent velocity 
components. The normal molecular cloud material at velocities 
of-13.5...-13kms-Mstracedby l^CO2-l,[CI],andCOll- 
10 in #1 and #3. In this complementary analysis approach, in- 
volving positions not all on the cut discussed above, we find that 
at #2, CO 11-10 is, significantly blueshifted, showing a local in- 
version of the velocity gradient similar to the inversion seen in 
[C ii] along the cut. This means that the embedded B 1 star not 
only creates a local H ii region, but it also aff'ects the dynamics 
of the hot gas in the PDR around the H ii region on a scale of 
about 30'', i.e. 0.1 pc. 

The [Cii] line is much wider at all positions. It includes the 
velocity of the molecular cloud but also an additional blueshifted 
component with velocities of about -16 knis~^ an d a wing ex- 



20 40 60 



tending down to -21 kms~ . Boreiko et al. (19901 showed that 
the -16 km s"^ component matches the vlsr of the Ha emission 
from the Hii region SI 55. At #1, this -16kms~^ is the domi- 
nant component, showing that most of the [C ii] emission stems 
from the H ii region there. The broad blue wing is, however, also 
present at #2 and #3, and actually strongest at #3. This may rep- 
resent ionized material that is ablating from the cloud surface. 
Depending on the exact geometry of cloud and main radiation 
source, we find two possible interpretations. If HD 217086, lo- 
cated slightly behind the cloud ( |Felli et al.|1978 ), represents the 
main sources of UV photons, the acceleration of low-density 
material through the radiation pressure could lead to the wing 
of blue shifted material. If HD 217061, in front of the cloud, is 
the main radiation source, we could be tracing the signat ure of 
a photo-evaporation flow from the surface as predicted bylGorti 
|& HoUenbach ( 2002| ). Detailed analysis confirming either of the 
two possibilities is beyond the scope of this paper. 



6. Conclusion 

The [C ii] emission is aff'ected on large scales by radiation from 
the Cepheus OB 3 association with a gradual change from ma- 
terial stemming mostly from the S155 Hii region in the NW 
through material photoevaporating from the molecular cloud 
surface, and since it is potentially accelerated through radiation 
pressure, to PDR material within the molecular cloud that is ion- 
ized through the FUV radiation, but resides at the velocities of 
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Fig. 5. (a) Map of the vlsr values (in kms"^) for Gaussian fits 
to the [C ii] spectra, (b) PV diagram for [C ii] along the cut indi- 
cated in (a). The white line shows the position of the [C ii] peak 
and the solid black line the centroid velocity along the cut. 
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Fig. 6. Spectra averaged over an area of 33'' diameter around the 
CO 11-10 peaks #1, #2, and #3. 



the parental cloud with the known large-scale velocity gradient. 
For the tracers of the denser material, we find no large-scale ef- 
fects but local velocity shifts due to the embedded ultracompact 
Hii regions. They create a partially inverted velocity gradient, a 
local broadening of the line profiles and disperses the material 
around them, but do not create any long-range eff'ects. 

Overall, we found no indications that a radiative trigger- 
ing mechanism creates the sequential star formation in Cepheus 
B. Radiative star- formation triggering may work on very small 
scales around the H ii regions mostly unresolved here, but cannot 
explain a large-scale sequence. 
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